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ABSTRACT 

We present the first measurement of the globular cluster populations of galaxies in Hickson 
compact groups, in order to investigate the effect of these high density environments on the 
formation and evolution of globular cluster systems. Based on V and R band images that we 
obtained of HCG 22a and HCG 90c with the ESO New Technology Telescope (NTT), we find a 
total globular cluster population of 1590 ± 854 for HCG 22a and 2136 ± 718 for 90c. The specific 
frequency for HCG 22a was found to be S N = 1.9 ± 1.0 and S N = 3.4 ± 1.1 for HCG 90c. A 
power-law fit to the globular cluster radial profile of HCG 22a yields a ~ /j-2.oi±o.30 anf j j or 
HCG 90c we found a ~ ^-i-20±o.i6 j^ com parison of the globular cluster radial profiles with 
the surface brightness of the parent galaxy shows that the globular cluster systems are at least 
as extended as the halo light. The measured values for the specific frequency are consistent with 
a scenario in which the host galaxies were in a low density "field-like" environment when they 
formed their globular cluster systems. 

Subject headings: galaxies: compact groups — galaxies: individual (NGC 1199, NGC 7173) — globular 
clusters: general 



1. Introduction 

Tremendous progress has been made in re- 
cent years in the detection of globular cluster 
(GC) populations around galaxies beyond the Lo- 
cal Group (e.g., Harris 1991; Blakeslee 1999). The 
globular cluster populations of well over 100 galax- 
ies up to distances > 100/i -1 Mpc have now been 



1 Based on observations collected at the European 
Southern Observatory, Chile. 



studied in various environments in the hope of 
gaining further insight into galaxy formation and 
evolution. Studies have shown that the number 
of globular clusters associated with any particu- 
lar galaxy varies widely, from none for M32 to 
- 20,000 for the cD galaxy M87 (e.g., Harris 
1991). 

A useful parameter for quantifying the total 
number of globular clusters per galaxy is the spe- 
cific frequency, S N = Ar tot 10 - 4 ( Mv+15 ), which is 



defined as the total number of globular clusters, 
Ntot, per unit galaxy luminosity (normalized to 
My — —15; Harris & van den Bcrgh 1981). Previ- 
ous studies have shown that globular cluster pop- 
ulations vary systematically with galaxy type and 
environment (Harris 1991; West 1994). Charac- 
teristic values of Sn range from Sjv < 1 for field 
spiral galaxies, to Sn — 2 — 4 for field ellipticals, 
to Sn — 5 — 7 for ellipticals in dense environments, 
and up to Sn — 10 — 20 for supergiant ellipticals at 
the centers of rich clusters. (We note that normal- 
ization with respect to bulge luminosity for spiral 
galaxies would increase the value of Sn, see Cote 
et al. 2000.) Yet, which factor has the greatest in- 
fluence on globular cluster formation — galaxy type 
or environment — is unclear. 

Several competing theories have been proposed 
to explain the origin of the observed variations 
in globular cluster systems of different galaxies. 
Among these are: galaxy mergers (Schweizer 1986; 
Ashman & Zepf 1992), accretion via tidal stripping 
or merging (Cote, Marzke, & West 1998; Cote et 
al. 2000), and two-stage collapse models (Forbes, 
Brodie, & Grillmair 1997). At present, it is un- 
clear whether a single mechanism, or a combina- 
tion of proposed formation scenarios, can explain 
the properties of globular cluster systems (Kissler- 
Patig 2000) . Determining which scenario most ac- 
curately describes the formation of globular clus- 
ters will require the investigation of GC systems 
in a wide range of environments. 

In this paper we present observations of the 
globular cluster population of two galaxies in Hick- 
son compact groups; HCG 22a (NGC 1199) and 
HCG 90c (NGC 7173). Compact groups are in- 
teresting because, although they are small (typ- 
ically 3-7 member galaxies per group), they are 
extremely dense systems, comparable in density 
to the cores of rich Abell clusters (e.g., Hickson 
1982; Hickson ct al. 1992). This makes them a 
useful tool for probing the influence of environ- 
ment on globular cluster formation. 

Figure 1 and 2 show V band images of HCG 22a 
and 90c, respectively. Both HCG 22 and 90 have 
been detected by ROSAT and each group has been 
shown to be a weak X-ray source (Ponman et al. 
1996). Table 1 summarizes the general properties 
of these systems. We discuss each in detail below. 



1.1. Hickson Compact Group 22 

HCG 22 initially was thought to contain five 
members (three ellipticals and two spirals) but 
later redshift measurements determined that two 
of the elliptical galaxies (HCG 22d and 22e) were 
most likely background galaxies since their helio- 
centric radial velocities are <~ 7000 km s _1 greater 
than the remaining group members (Hickson et al. 
1992). The median heliocentric radial velocity of 
the three bona fide group members is 2700 km s , 
with a velocity dispersion of 43.7 km s . 

From Table 1, the heliocentric velocity of HCG 
22a is 2705 km s _1 which indicates that this galaxy 
is most likely located near the group center of 
mass. Hickson et al. (1992) estimated the mass-to- 
light ratio of HCG 22 to be M/L ~ I.S/T^Mq/Lq 
and a projected median galaxy-to-galaxy separa- 
tion of 26.9/i _1 kpc. The group crossing time, de- 
fined as the ratio of the crossing time to the age 
of the universe, is Ht c = 0.1905 which is greater 
than the median value of Ht c — 0.016 determined 
for all of the Hickson compact groups combined 
(Hickson et al. 1992). This suggests that galaxy 
mergers and interactions may be less important 
for HCG 22 than for other compact groups with 
shorter crossing times. 

1.2. Hickson Compact Group 90 

HCG 90 consists of four bright galaxies (two el- 
lipticals, a spiral, and an irregular) with a median 
heliocentric group radial velocity of ~ 2640 km s _1 
and a velocity dispersion of ~ 100 km s _1 (Hick- 
son et al. 1992). The mass-to-light ratio of HCG 
90 is given by Hickson et al. (1992) as M/L ~ 
12.3h~ 1 M Q /L ( 7 ) and a median projected galaxy 
separation of 29.5/i~ 1 kpc. Table 1 lists the helio- 
centric radial velocity of HCG 90c as 2696 km s _1 , 
which indicates that this galaxy is probably close 
to the center of mass of the group. 

The group crossing time was estimated by Hick- 
son et al. (1992) to be Ht c = 0.0224 , which is 
only ~ 12% of the estimated value for HCG 22. 
This suggests that galaxy mergers and interactions 
are likely to have been more important in HCG 
90 than for HCG 22 and it is interesting to note 
that HCG 90 contains a galaxy (HCG 90d) that 
is clearly interacting with other group members 
(Longo et al. 1994; Plana et al. 1999). 

In this study, we have measured the globular 



cluster specific frequency of HCG 22a and HCG 
90c. In § 2 we present the observations and data 
reductions. Section 3 presents the results for HCG 
22a and 90c in terms of their globular cluster lumi- 
nosity functions, radial distributions, and specific 
frequencies. In § 4 we discuss the results of this 
study in the context of globular cluster formation 
and the nature of Hickson compact groups. 

Throughout this paper, unless otherwise speci- 
fied, we use Hq = 50 km s _1 Mpc _ . 

2. Observations and Data Reductions 

Observations of HCG 22a and 90c were ob- 
tained on the nights of October 9 and 10, 1993 
using the 3.5m New Technology Telescope (NTT) 
operated by the European Southern Observatory 
(ESO) in La Silla, Chile. All images were taken 
with the ESO Multi-Mode Instrument (EMMI) 
with a Loral 2048 CCD mounted on the red arm. 
The Loral 2048 CCD chip consists of 2048x2048 
15/im pixels, giving a scale of 0.29" per pixel at the 
f/11 Nasmyth focus. Technical limitations during 
the observing run limited the effective area of the 
chip to 1700x1700 pixels, resulting in a projected 
sky coverage of 8.2' x 8.2'. 

The observations consisted of 5 x 900 second 
exposures in the Cousins V and R passband for 
HCG 22a, and 6 x 900 seconds in V and R for 
HCG 90c. Since globular clusters around these 
galaxies can be detected only as a statistical ex- 
cess of star-like objects, a series of 5 x 900 second 
exposures in V and R were also taken of a control 
field located within 3° of the target galaxies in or- 
der to ascertain the expected number of contami- 
nating Galactic stars and unresolved faint galaxies 
along the line of sight. The seeing varied during 
the observing run from 0.87" to 0.98" (FWHM), 
as measured from the individual galaxy frames. 

Processing of the raw images was conducted 
within the IRAF 2 environment and involved the 
standard procedure of bias correction and flat 
fielding using dome flats. After this initial pro- 
cessing, each set of individual galaxy and control 
field images were registered to their common co- 
ordinate system and median combined to yield 



2 IRAF is distributed by National Optical Astronomy Obser- 
vatories, which arc operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



a higher S/N image. Instrumental magnitudes 
were transformed to the standard UBVRI system 
by measuring several standard stars from Landolt 
(1992). 

At the distance of the target galaxies (~ 
54 Mpc), and with the typical seeing of the com- 
bined frames (~ 0.9"), individual globular clusters 
are unresolved and appear as an apparent excess 
of star-like objects around their parent galaxy. In 
order to maximize the detection of individual glob- 
ular clusters and to perform accurate photometry, 
the light from each host galaxy was removed. This 
was done by using the STSDAS package in IRAF 
to model the light distribution of each galaxy using 
the tasks ellipse and bmodel, which fit isophotcs 
and then generate models of each galaxy's surface 
brightness profile. The constructed galaxy models 
were then subtracted from the appropriate image. 

The galaxy-subtracted image of HCG 22a re- 
vealed a large, central, diffuse structure which is 
the residual of a known central dust lane (Franx, 
Illingworth, & Heckman 1989). In order to remove 
this feature, the galaxy-subtracted image was me- 
dian filtered using a ring median filter (Seeker 
1995) having an inner radius of 6.36 pixels and 
a width of 1 pixel. The resultant smoothed image 
was subtracted from the previous image, which 
helped to remove most of the dust lane feature. 
A second rectangular median filter of dimension 
10x10 pixels was applied to the resultant image 
and this smoothed image was subtracted from the 
previous image. The galaxy light removal proce- 
dure was applied equally to the V and R images. 

The HCG 90c galaxy model was subtracted 
from the parent image and the result was median 
filtered using a ring filter with the same dimension 
as that used for HCG 22a. The resultant smoothed 
image was then subtracted from the previous im- 
age, adequately removing all galaxy light. This 
procedure was applied equally to the V and R im- 
ages. 

Object detection and classification was con- 
ducted using the Faint Object Classification and 
Analysis System (FOCAS) within the IRAF envi- 
ronment (Jarvis & Tyson 1981; Valdes 1993). A 
threshold of 2.5cr above background and a min- 
imum object area of six pixels was used for the 
detection criteria. The detected objects were clas- 
sified by FOCAS as star-like or galaxy-like by com- 
paring a model PSF, generated by the average 



profile of a minimum of ten isolated, unsaturated 
stars, with the light profile of each object. The ob- 
ject list was then culled of saturated and galaxy- 
like objects. A master catalog of stellar-like ob- 
jects detected in both the V and R images was gen- 
erated for each galaxy and the control field. Only 
objects that were matched in both filters were in- 
cluded in the final sample, thereby greatly reduc- 
ing the possibility of false detections of objects. 
DAOPHOT was used to perform the photometry 
on all detected objects. 

Extinction corrections (Burstein & Heiles 1982) 
were determined to be negligible for HCG 90 and 
the control field (A B = 0.0). For HCG 22, we 
determined Ay = 0.08 magnitude and Ar = 0.06 
magnitude. 

2.1. Completeness Function 

Since the globular cluster population of HCG 
22a and 90c are detected only at very faint mag- 
nitudes (V > 22), the probability of detecting an 
object as a function of magnitude must be deter- 
mined as accurately as possible. A "complete- 
ness" function was calculated by adding artifi- 
cial stars (generated using the DAOPHOT routine 
Addstar) to the galaxy and control field images 
and then re-running the FOCAS routines. The 
fraction of detected artificial stars as a function of 
magnitude was then used to derive the complete- 
ness function. 

Artificial stars were added to both the V and 
R images using the same spatial coordinates and 
with V — R = 0.4 (the approximate color of Milky 
Way globular clusters; see Peterson 1993). FO- 
CAS was then run exactly in the same manner as 
that used for real objects, including keeping only 
those objects detected in both the V and R images. 

The photometric magnitude limit of the data is 
defined as the magnitude at which the probability 
for detecting an object equals 50% (Harris 1990) 
. For magnitudes fainter than this limit, object 
information is much less reliable and thus all such 
objects were discarded from further analysis. The 
completeness limit for HCG 22a was determined 
to be V = 24.9. 

Since noise and degree of crowding increases 
towards the center of the subtracted galaxy, the 
completeness function will vary with radial dis- 
tance from the galaxy center. Artificial stars were 



placed at different radii from the galaxy center and 
the completeness limit, as a function of radial dis- 
tance, was measured. The completeness function 
was measured as a function of both magnitude and 
position from the center of the subtracted galaxy. 

For HCG 22a, it was determined that for radii < 
70 pixels, the completeness limit was several mag- 
nitudes brighter than for radii > 70 pixels (part of 
this is due to the small disk-like artifact found at 
the center of HCG 22a). All objects found within 
70 pixels of the galaxy center were discarded and 
this inner region was not used in any further anal- 
ysis. For the area outside of this region, the com- 
pleteness function was found to be approximately 
the same. 

For HCG 90c, the area within 60 pixels of the 
galaxy center was determined to be unusable due 
to its bright magnitude completeness limit. This 
region was masked out and not used in any further 
analysis. The completeness limit of the area out- 
side of this region was determined to be V — 24.8. 

For the control field, artificial star experiments 
indicate a completeness limit of V = 24.6. 

Photometric errors were calculated by deter- 
mining the standard deviation of the magnitude 
difference between added and recovered artificial 
stars. The root-mean-square of the random er- 
ror for HCG 22a ranged from a rms s=s 0.01 at 
V = 22 to <r rms w 0.06 at V = 24.6. For HCG 
90c, a rms w 0.01 at V = 22 and increased to 
a rms w 0.09 at V = 24.6. 

2.2. Selection Criteria 

The detection of globular cluster candidates can 
be improved using criteria such as magnitude and 
shape parameters which are characteristic of glob- 
ular clusters in general. Since objects detected 
in both galaxy fields must be compared to ob- 
jects detected in the control field, all objects in 
the galaxy image catalogs with magnitudes fainter 
than V = 24.6 (the completeness limit of the con- 
trol field) were culled. Also, at the distance of 
HCG 22a and 90c, the number of globular clusters 
brighter than V = 22 is negligible and hence all 
objects brighter than this limit were removed from 
the object catalogs. 

Besides the magnitude criteria, all objects 
that lie within the masked regions (such as near 
the galaxy centers and regions containing bright 



galaxies) were eliminated from the object catalogs. 
In order to improve the S/N of detected globular 
clusters over background galaxies, all objects de- 
tected more than 2.13' from the center of HCG 
22a were discarded (at this radius the density of 
globular cluster candidates becomes negligible; see 
section 3.2). For HCG 90c, a cutoff radius of 2.58' 
was chosen. 

Figures 3 and 4 depict the spatial distribution 
of globular cluster candidates for HCG 22a and 
HCG 90c. A concentration of star-like objects sur- 
rounding each galaxy center is clearly seen. 

Figures 5, 6, and 7 give the color-magnitude 
distribution for the galaxy and control fields. It 
appears from these figures that the distribution 
of stellar-like objects in HCG 22a and HCG 90c 
are generally more concentrated at a specific color 
than those objects found in the control field. Due 
to the small number of detected objects and the 
very short color baseline, we have not imposed a 
color selection criteria. 

3. Results 

3.1. Globular Cluster Luminosity Func- 
tion 

At a limiting magnitude of V ~ 24.6 we sam- 
ple only the brightest fraction of globular clusters 
present around HCG 22a and 90c. Consequently, 
it is necessary to extrapolate counts to include the 
more numerous faint globulars not detected in our 
images. 

Globular cluster luminosity functions are uni- 
versally well described by a Gaussian function 
(Harris 1991) of the form 



N(V) = A exp 



— (m — m o y 



2ct 2 



(1) 



where N(V) is the number of globular clusters, A 
is the amplitude, m is the apparent magnitude of 
individual globular clusters, m is the magnitude 
of the turnover, and a is the dispersion. 

The globular cluster luminosity function (GCLF) 
for each galaxy was constructed from the statisti- 
cal subtraction of star-like objects measured from 
the control field. Unresolved objects in the galaxy 
field will consist of true globular clusters, fore- 
ground stars, and faint background galaxies which 
have stellar-like profiles and were miss-classified 



by FOCAS. In fact, most of the contamination 
will be due to faint unresolved galaxies since one 
expects less than one foreground star per square 
arcminute at V > 21 (e.g., Ratnatunga & Bahcall 
1985). Since object detection and classification 
were carried out identically for both the galaxy 
and control fields, the difference in the number of 
star-like objects between the galaxy fields and the 
control field (corrected for completeness and nor- 
malized with respect to spatial coverage) should 
be, statistically, globular cluster candidates. 

Fitting a Gaussian function to the GCLF data 
involves simultaneously fitting three parameters; 
the turnover, dispersion, and the amplitude. How- 
ever, since the GCLF data for HCG 22a and 90c 
are highly uncertain near the completeness limit, 
it was decided not to attempt to simultaneously 
solve for all three parameters. The absolute mag- 
nitude of the turnover of the GCLF is observed 
to be roughly the same for all galaxies, regard- 
less of morphological type or local environment 
(Harris 1991; Jacoby et al. 1992). We have there- 
fore adopted My = —7.4 ± 0.2 for the absolute 
magnitude of the turnover, which is consistent 
with previous studies (e.g., Fleming et al. 1995). 
Since the dispersion of the Gaussian function has 
been shown to have only a narrow range of val- 
ues centered on a gauss = 1.4 (Harris 1991; Ja- 
coby et al. 1992), we have decided to adopt this 
value. A non-linear weighted least-squares fit of 
the GCLF of HCG 22a, using M v = -7.4 ± 0.2, 
&gauss = 1-4 and a distance of 54.1 Mpc, gives 
A = 38.4 ± 8.1 with xl = 1-4. For HCG 90c, we 
found A = 51.3 ± 7.6 and xt = 1 -9, using a dis- 
tance of 53.9 Mpc. Figure 8 and 9 show the GCLF, 
with the best-fit Gaussian function, for HCG 22a 
and HCG 90c, respectively. Counts of stellar-like 
objects in the galaxy and control fields are given 
in Table 2 and 3 for HCG 22a and HCG 90c, re- 
spectively. 

3.2. Radial Distribution 

In addition to correcting globular cluster counts 
for incomplete sampling of the GCLF as discussed 
above, it is also necessary to correct the counts for 
the limited areal coverage around each galaxy. 

The globular cluster radial distribution was 
measured by dividing the galaxy images into sev- 
eral annuli centered on each galaxy. For each an- 
nulus, the total number of globular cluster candi- 



dates were corrected for incompleteness and sub- 
tracted by the properly normalized control field 
counts. 

The radial distribution was characterized using 
a power-law of the form, a = AR a (e.g., Flem- 
ing et al. 1995). Using weighted non- linear least- 
squares, the best-fit values for the radial distri- 
bution of HCG 22a was determined to be A = 
7.2 ± 1.28 and a = -2.01 ± 0.30, with xl = 1-3. 
Figure 10 depicts the radial distribution of globu- 
lar cluster candidates for HCG 22a along with the 
best-fit power-law. 

For HCG 90c, the best-fit values for the radial 
distribution was found to be A = 13.99 ± 1.42 
and a = -1.20 ± 0.16, with xl = 1-5. Figure 11 
depicts the radial profile of the globular cluster 
candidates, along with the best-fit power-law. 

A useful comparison can be made between the 
spatial extent of the globular cluster system and 
the stellar halo of the parent galaxy. Figure 12 
shows the comparison of the radial profile of the 
globular cluster system and the surface brightness 
for HCG 22a. The surface brightness data has 
been taken from Franx et al. (1989) and has been 
scaled vertically to match the GC profile near its 
center. A fit to the radial dependence of the halo 
surface brightness, /x, gives, \i ~ R~ 23 , which is 
consistent with the slope of the radial profile of the 
GCs (a = -2.01 ±0.3). Thus, the globular cluster 
system is as spatially extended as the halo light of 
HCG 22a, which is consistent with other studies 
where the GC profiles are either similar to or more 
extended than the galaxy halo (e.g., Fleming et al. 
1995). 

The surface brightness profile of HCG 90c was 
determined from Penereiro (1994) and is compared 
with the radial distribution of its globular clus- 
ter population in Figure 13. The surface bright- 
ness profile has been scaled vertically to match 
the globular cluster distribution. A fit to the halo 
profile gives a slope of /i ~ R~ 15 and hence, the 
globular cluster population is spatially more ex- 
tended than the galaxy halo. We also note that 
the halo light of HCG 90c is not well defined due 
to the high surface brightness of the intra-group 
light (White et al. 2001). 



3.3. Specific Frequency 

As mentioned in Section 1, a useful way of quan- 
tifying the globular cluster population of different 
galaxies is to compare values of specific frequency, 
S N = N tot W Ai * Mv+15 \ The first step in calculat- 
ing Sn is to determine the total number of globu- 
lar clusters, N to t- The power- law radial distribu- 
tion function, a = AR a , must be integrated from 
an inner to an outer radius in order to account for 
the globular clusters that were missed due to spa- 
tial incompleteness. Also, the number of globular 
clusters must be corrected for incomplete coverage 
of the luminosity function. 

3.3.1. Hickson Compact Group 22a 

In order to integrate the power-law radial dis- 
tribution function, inner and outer radial limits 
must be chosen. Globular cluster distributions 
normally do not extend to radii inward of 1 kpc 
from the center of the parent galaxy (e.g., Harris 
1991; Bridges, Hanes, & Harris 1991), the most 
likely cause being tidal shocking and dynamical 
friction which have destroyed globular clusters at 
distances close to the galaxy center (e.g., Wein- 
berg 1993). For an outer radius, globular cluster 
systems generally do not extend past 50 kpc from 
the center of a typical elliptical galaxy (e.g., Har- 
ris & van den Bergh 1981; Hanes & Harris 1986; 
Fleming et al. 1995). Thus for this study, an inner 
radius of 1 kpc and an outer radius of 50 kpc was 
adopted for HCG 22a and 90c. These radial lim- 
its for HCG 22a correspond to angular distances 
of 0.06' and 3.1'. A consistency check is provided 
by the fact that the radial counts for HCG 22a 
and 90c appear to reach background levels approx- 
imately 50 kpc from the galaxy center. 

The integral of the power-law density profile 
gives N to t — J o~ c i 2irR dR, where a c i represents 
the number density of globular clusters. For HCG 

22a we have N tot = 2ttA j^ 8 R a+1 dR. Solving 
this equation gives the total number of globular 
clusters, N tot — 179 ± 54 before including the cor- 
rection for faint globular clusters that are below 
our detection limit. If the outer radius was set to 
a smaller value, say 25 kpc, the total number of 
globular clusters would be ~ 148 and hence, N to t 
is not extremely sensitive to our chosen value of an 
outer radius for a steep radial slope of a ~ —2.0. 

After correcting for the radial incompleteness of 



detected globular clusters, the incomplete cover- 
age of the luminosity function must be accounted 
for. Since the luminosity function was assumed to 
be a Gaussian, the area under the Gaussian func- 
tion that was accounted for will give a measure of 
the percentage of observed globular clusters. In- 
tegrating the curve from — oo to the photometric 
limit of the data (V — 24.6), indicates that 11±3% 
of the distribution was observed. Correcting N to t 
for the unobserved part of the luminosity function 
gives N t c ° t r = 1590 ± 854. 

Once the total number of globular clusters was 
determined, the specific frequency was calculated 
from the absolute V magnitude of the parent 
galaxy. Using an absolute magnitude of My = 
-22.31 ± 0.07 as given by NED 3 , the specific fre- 
quency was determined to be Sn = 1-9 ± 1.0 for 
HCG 22a (using H Q = 75 km s^Mpc -1 , we find 
S N = 1.7 ±0.6). 

3.3.2. Hickson Compact Group 90c 

The specific frequency for HCG 90c was calcu- 
lated in the analogous way as for HCG 22a. The 
total number of globular clusters, for the chosen 
radial limits of 1 kpc to 50 kpc (0.07' to 3.5'), was 
determined to be N to t = 270 ± 28. Changing the 
outer radius to 25 kpc yields N ta t ~ 150 and thus 
the total number of globular clusters is sensitive 
to the choice of an outer radius because of the 
shallow radial profile with a ~ —1.2. 

Correcting for the unobserved part of the lu- 
minosity function, we determined that 13 ± 3% 
of the total number of globular clusters were ob- 
served and that N™£ = 2136 ± 718. Using an 
absolute magnitude of My = —21.99 ± 0.03, we 
find that S N = 3.4 ± 1.1 for HCG 90c (using 
Ho = 75 km s _1 Mpc _1 , we find S N = 4.4 ± 1.3). 

4. Discussion 

In this paper we have investigated the globu- 
lar cluster population of two prominent elliptical 
galaxies, each located near the center of a Hick- 
son compact group. The globular cluster popula- 
tion has been characterized by measuring the num- 
ber of globular clusters per unit galaxy luminosity 
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ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. 



(Sn), radial distributions, and the globular cluster 
luminosity function. Although our sample is lim- 
ited to only two compact group galaxies, we can 
draw some tentative conclusions. 

The most important aspect of this study is the 
determination of the specific frequency of HCG 
22a and HCG 90c. As mentioned in section 1, 
specific frequency values for field ellipticals gener- 
ally range from Sn = 2 — 3, in contrast to values 
measured for cD galaxies at the center of rich clus- 
ters (Sn = 10— 20). Although Sn seems to follow 
a general trend that as the local density increases 
the number of globular clusters per unit galaxy lu- 
minosity also increases, there are exceptions (e.g., 
Harris, Pritchet, & McClure 1995; Woodworth & 
Harris 2000). The evidence seems to suggest that 
a high density environment, at least at the current 
epoch, may be a necessary but not a sufficient con- 
dition for the formation of high Sn galaxies. The 
relatively low values of Sn for HCG 22a and 90c 
(S N = 1.9 ± 1.0 and S N = 3.4 ± 1.1, respectively) 
suggests that these galaxies formed their globu- 
lar clusters at a time when they were located in a 
"field- like" low density environment. 

The exact nature of Hickson compact groups 
has been controversial since some studies have sug- 
gested that the majority of compact groups are 
chance projections of large filamentary structures 
(Hcrnquist, Katz, & Weinberg 1995) or a superpo- 
sition of galaxies in a much larger group or poor 
cluster (Walke & Mamon 1989). N-body simula- 
tions have indicated that, given the small cross- 
ing times and low velocity dispersions (Hickson et 
al. 1992), group member galaxies would combine 
to form a single elliptical galaxy on the order of 
a few crossing times (e.g., Barnes 1985; Mamon 
1987; but see Mamon 2000) and several examples 
of isolated ellipticals which may be the surviv- 
ing remnants of galaxy groups have been discov- 
ered (Mulchaey & Zabludoff 1999). This suggests 
that if compact groups are gravitationally bound, 
group members may be just beginning to come to- 
gether at the present epoch and that during the 
time of globular cluster formation, these galaxies 
were situated in a lower density environment. This 
scenario is consistent with the specific frequency 
values measured for HCG 22a and 90c. 

Evidence to support the hypothesis that these 
groups are bound structures include; the detec- 
tion of Hj and X-rays (Williams & Rood 1987; 



Pildis, Brcgman, & Evrard 1995; Ponman et al. 
1996), indicating significant dynamical evolution 
and galaxy interaction, the decrease in the rela- 
tive number of spirals with decreasing group cross- 
ing times (Hickson et al. 1992), and the presence 
of tidal tails and debris in some compact groups 
(Hunsberger, Charlton, & Zaritsky 1996). 

As mentioned in section 1, both HCG 22 and 
90 arc weak X-ray emitters. Radial velocity mea- 
surements imply that both HCG 22a and 90c are 
located near the bottom of their group potential 
well. Also, HCG 90b is clearly interacting with 
HCG 90d (Plana et al. 1999). HCG 90 has also 
been shown to contain a spatially extended, dif- 
fuse, light component (White et al. 2001). These 
facts suggests that both HCG 22 and 90 are grav- 
itationally bound structures. 

In recent years, it has been shown that the glob- 
ular cluster systems of many large elliptical galax- 
ies contain two or more chemically distinct sub- 
populations, usually inferred on the basis of color 
information (e.g., Ostrov, Geisler, & Forte 1993; 
Geisler, Lee, & Kim 1996; Cote et al. 1998). A 
number of different theories have been proposed 
to explain the origin of these different globular 
cluster populations; among these are models in 
which two or more bursts of globular cluster for- 
mation occur (e.g., Ashman & Zcpf 1992; Forbes 
& Forte 2000) or accretion of globular clusters 
from other galaxies (e.g., Cote et al. 1998, 2000). 
Given the rather unique environments represented 
by compact groups, and the roles that mergers 
may have played, it would be interesting to as- 
certain whether Hickson compact group galaxies 
also exhibit multiple globular cluster populations. 
Unfortunately, the V-R information provided by 
this investigation does not give us a sufficiently 
long color baseline to explore this possibility. 

5. Conclusions 

We have studied the globular cluster systems 
around the early-type galaxies HCG 22a (NGC 
1199) and 90c (NGC 7173) using V and R band 
data obtained from ESO's NTT. Globular clus- 
ters have been detected as a statistical excess of 
star-like objects surrounding each galaxy. Object 
detection and classification was performed using 
FOCAS on galaxy-subtracted images. 

The specific frequency of each galaxy was mea- 



sured assuming a turnover and dispersion of the 
globular cluster luminosity function. We found 
that for HCG 22a, S N = 1.9 ± 1.0, and for HCG 
90c, S N = 3.4 ± 1.1. The total number of glob- 
ular clusters for HCG 22a was determined to be 
N tat = 1590 ± 854 and for HCG 90c we found 
Ntot = 2136 ± 718. We also derived the globular 
cluster surface density profiles and found that for 
HCG 22a, the underlying starlight is as extended 
as the GC system. For HCG 90c, we found that 
the GC system is more extended than the galaxy 
halo. 

From the results of this study, it is clear that 
further observations are needed to help under- 
stand the formation of globular clusters and the 
nature of Hickson compact groups. Observations 
which clearly reach the turnover in the globular 
cluster luminosity function would be of great value 
since they would provide a more robust determina- 
tion of the GCLF, radial distribution, and specific 
frequency. Observations design to detect any bi- 
modality in the color distribution would provide 
invaluable information on the past history of the 
globular cluster systems. 

Additional observations of a larger sample of 
compact group galaxies would allow the investiga- 
tion of potential correlations between galaxy lumi- 
nosity and Sn (since the local environment would 
be the same for a pair of galaxies in the same 
group) or other factors such as X-ray luminosity, 
mass-to-light ratio, and distance from the group 
dynamical center. 
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Fig. 1. — V band image of Hickson compact group 
22a (NGC 1199). 

Fig. 2. — V band image of Hickson compact group 
90 with HCG 90c located near the center position. 



Fig. 3. — Positions of globular cluster candidates 
surrounding HCG 22a with V < 24.6. The central 
70 pixel masked-out region is shown along with 
the galaxy center at X = 885 and Y = 647 pixels. 
Only objects having radii 70 < R < 441 pixels are 
shown. Note the concentration of star-like objects 
around the center of HCG 22a. 

Fig. 4. — Positions of globular cluster candidates 
surrounding HCG 90c with V < 24.6. The central 
60 pixel masked-out region is shown along with 
the galaxy center at X = 887 and Y = 652 pixels. 
Only those objects having radii 60 < R < 535 pix- 
els are shown, excluding the area in the lower left 
quadrant which was masked out due to the pres- 
ence of HCG 90b and d. Note the concentration 
of star-like objects around the center of HCG 90c. 

Fig. 5.— Color-magnitude diagram of the globular 
cluster candidates surrounding HCG 22a. Note 
the enhanced concentration of stellar-like objects 
near V — R = 0.6. 

Fig. 6. — Color- magnitude diagram of the globular 
cluster candidates surrounding HCG 90c. Note 
the enhanced concentration of stellar-like objects 
near V — R = 0.5. 

Fig. 7. — Color-magnitude diagram of stellar-like 
objects in the control field. Note that objects ap- 
pear to be scattered to a greater extent than those 
found for HCG 22a and 90c. 
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Fig. 8. — Globular cluster luminosity function of 
HCG 22a. The best-fit Gaussian function is shown 
as the dashed line with a turnover of my — 26.3 ± 
0.2 and a dispersion of a = 1.4. 

Fig. 9. — Globular cluster luminosity function of 
HCG 90c. The best-fit Gaussian function is shown 
as the dashed line with a turnover of my = 26.0 ± 
0.2 and a dispersion of a = 1.4. 

Fig. 10. — Radial distribution of globular clusters 
around HCG 22a. The best-fit power-law func- 
tion, a ~ ^-2.oi±o.30^ j g s j 10wn as the dashed line. 
The radius is given as the geometric mean of the 
individual annulii,< R > = y ' Ri n Rout, and a is 
the number of star-like objects per arcminute . 

Fig. 11. — Radial profile of globular clusters sur- 
rounding HCG 90c. The best-fit power-law, a <~ 
jj-i.20io.i6 i s shown as a dashed line. The radius 
and a are as defined for HCG 22a. 

Fig. 12. — The radial profile of the globular cluster 
system of HCG 22a (dashed line) compared to the 
R-band galaxy surface brightness from Franx et 
al. (1989) (solid line). The surface brightness pro- 
file has been scaled vertically in order to match 
the GC profile near the center. From the figure, 
it is clear that the globular cluster system is as 
spatially extended as the halo light of the parent 
galaxy. The slope of the power-law fit to the GC 
profile is a = — 2.01 ± 0.3 and for the halo light 
the slope is a ~ —2.3. 

Fig. 13. — The radial profile of the globular clus- 
ter system of HCG 90c (solid line) compared to the 
R-band galaxy surface brightness from Pcnereiro 
(1994) (dashed line). The surface brightness pro- 
file has been scaled vertically in order to match 
the GC profile near the center. From the figure, 
it is clear that the globular cluster system is more 
spatially extended than the halo light of the par- 
ent galaxy. The slope of the power-law fit to the 
GC profile is a = —1.20 ± 0.16 and for the halo 
light the slope is a ~ —1.5. 
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Table 1 
General properties of HCG 22a and HCG 90c 





HCG 22a 1 


HCG 90c 


a (1950) 


03 01 18.2 


21 59 08.8 


5 (1950) 


-15 48 30 


-32 12 58 


I (deg) 


199.22 


14.98 


b (deg) 


-57.31 


-53.08 


Hubble Type 


E2 


E0 


V (km s _1 ) 


2705 


2696 


S T (mag) 


12.24 


12.73 


B-R (mag) 


1.62 


2.23 



Tabic information from Hickson (1994). 
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Table 2 
Globular Cluster Counts: HCG 22a 



V 


N ga i. 


N bkg 


N 


22.1 


3.0 ± 1.7 


0.6 ±0.4 


2.4 ± 1.8 


22.3 


3.0 ± 1.7 


2.2 ±0.8 


0.8 ± 1.9 


22.5 


0.0 ±0.0 


0.8 ±0.5 


-0.8 ±0.5 


22.7 


5.0 ±2.2 


1.4 ±0.6 


3.6 ±2.3 


22.9 


7.0 ±2.6 


2.8 ±0.9 


4.2 ± 2.8 


23.1 


3.1 ± 1.8 


0.6 ±0.4 


2.5 ± 1.8 


23.3 


9.3 ±3.1 


1.6 ±0.3 


7.6 ±3.1 


23.5 


8.1 ±2.9 


2.8 ±0.9 


5.4 ±3.0 


23.7 


14.2 ±4.2 


4.4 ± 1.1 


9.7 ±4.3 


23.9 


16.1 ±4.4 


7.2 ± 1.4 


9.0 ±4.6 


24.1 


10.6 ±3.4 


10.0 ± 1.6 


0.6 ±3.8 


24.3 


35.4 ±6.6 


8.8 ± 1.7 


26.6 ±6.8 


24.5 


30.2 ± 7.0 


12.0 ± 2.0 


18.2 ±7.3 



Note. — Column 1 gives the mid-bin V mag- 
nitude. Columns 2 and 3 gives the raw counts 
of stellar-like objects in the galaxy and control 
field. Column 4 gives the background-corrected 

counts of stellar-like objects in the galaxy field. 
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Table 3 
Globular Cluster Counts: HCG 90c 



V 


N ga i_ 


Nbkg 


N 


22.1 


6.0 ±2.4 


0.6 ±0.4 


5.4 ±2.5 


22.3 


6.0 ±2.4 


2.5 ±0.9 


3.5 ±2.6 


22.5 


8.0 ±2.8 


0.9 ±0.5 


7.1 ±2.9 


22.7 


6.0 ±2.4 


1.6 ±0.4 


4.4 ± 2.5 


22.9 


13.0 ±3.6 


3.1 ± 1.0 


9.9 ±3.7 


23.1 


15.0 ±3.9 


0.6 ±0.4 


14.4 ±3.9 


23.3 


6.0 ±2.4 


1.9 ±0.8 


4.1 ± 2.6 


23.5 


10.0 ±3.2 


3.1 ± 1.0 


6.9 ±3.3 


23.7 


20.0 ±4.5 


5.0 ± 1.2 


15.0 ±4.6 


23.9 


22.0 ±4.7 


8.1 ± 1.6 


13.9 ±5.0 


24.1 


26.0 ±5.1 


11.2 ± 1.9 


14.8 ±5.4 


24.3 


50.1 ±7.4 


9.9 ± 1.9 


40.2 ± 7.7 


24.5 


38.8 ±6.9 


13.5 ±2.3 


25.4 ± 7.2 



Note. — Column 1 gives the mid-bin V mag- 
nitude. Columns 2 and 3 gives the raw counts 
of stellar-like objects in the galaxy and control 
field. Column 4 gives the background-corrected 
counts of stellar-like objects in the galaxy field. 
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